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ABSTRACT 


Modem  avionics  equipment,  such  as  super-resolution  direction-finding  systems, 
now  require  resolutions  on  the  order  of  20  to  22  bits.  Oversampled  analog-to-digital 
converter  architectures  offer  a  means  of  exchanging  resolution  in  time  for  that  in 
amplitude  and  represent  an  attractive  approach  to  implementing  precision  converters 
without  the  need  for  complex  precision  analog  circuits.  Using  oversampling  techniques 
based  on  sigma-delta  modulation,  a  convenient  tradeoff  exists  between  sampling  rate  and 
resolution.  One  of  the  major  advantages  of  integrated  optics  is  the  capability  to 
efficiently  couple  wideband  signals  into  the  optical  domain.  Typically,  sigma-delta 
processors  require  simple  and  relatively  low-precision  analog  components  and  thus  are 
well  suited  to  integrated  optical  implementations.  This  thesis  reviews  the  current  sigma- 
delta  methodology,  the  advantages  of  optical  integrated  circuits  and  presents  the  design  of 
a  second-order,  integrated  optical  sigma-delta  modulator.  Simulation  results  for  both  a 
first  and  second  order  architecture  are  presented  by  evaluating  the  transfer  characteristics 
numerically.  Design  parameters  such  as  limit  cycles  are  quantified  and  explained. 
Performance  issues  and  future  efforts  are  also  considered. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Analog-to-digital  converters  (ADCs)  are  basic  building  blocks  for  a  wide  variety 
of  digital  systems.  A  partial  list  of  ADC  applications  includes  process  control,  automatic 
test  equipment,  video  signal  acquisition,  audio  recordings  for  compact  disc  and  interfaces 
for  personal  computers.  There  exists  a  variety  of  approaches  to  the  ADC  design.  One 
approach,  known  as  delta  modulation,  involves  the  use  of  oversampling  methods.  First 
introduced  in  the  1 940s,  delta  modulation  uses  oversampling  and  single-bit  code  words  to 
represent  the  analog  signals  [Ref.  1].  The  simplest  approach  counted  the  output  bits  from 
the  delta  modulator  with  a  high  bit  representing  a  +1  and  a  low  bit  a  -1.  The  output  was 
then  resampled  at  the  Nyquist  rate.  Resolution  proved  to  be  a  problem,  since  achieving 
adequate  reproduction  of  speech  signals  required  oversampling  ratios  of  the  order  5,000. 
More  effective  digital  filtering  was  needed  to  prevent  the  high-frequency  modulation 
noise  from  aliasing  into  the  signal  band  when  it  was  resampled  at  the  Nyquist  rate. 

Unfortunately  at  that  time,  digital  filters  used  for  this  purpose  were  prohibitively 
expensive.  Candy  proposed  an  interpolative  technique  for  digital  filtering  [Ref.  2].  The 
idea  was  to  digitize  the  signal  through  the  use  of  a  coarse  quantizer  and  to  cause  the 
output  to  oscillate  between  the  quantized  levels  at  high  speed  so  that  its  average  value 
over  the  Nyquist  interval  was  an  accurate  representation  of  the  sampled  value.  The 
digital  filters  used  to  generate  this  average  were  inexpensive.  On  the  other  hand,  these 
digital  filters  also  proved  to  be  reliable  and  fairly  tolerant  of  circuit  imperfections.  The 
quantizers  for  these  interpolating  converters  utilized  a  noise-shaping  technique  which 
measures  the  quantization  error  in  one  sample  and  subtracts  it  from  the  next  input  sample 
value  [Ref.  3].  The  most  popular  form  of  this  noise-shaping  technique  is  known  as  sigma 
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delta  modulation.  Sigma  delta  modulators  employ  integration  and  feedback  in  iterative 
loops  to  obtain  high-resolution  A/D  conversions. 

Specifically,  a  sigma  delta  modulator  (IAM)  consists  of  an  analog  filter  and  a 
quantizer  enclosed  in  a  feedback  loop  [Ref.  4],  Together  with  the  filter,  the  feedback 
loop  acts  to  attenuate  the  quantization  noise  at  low  frequencies  while  amplifying  the  high- 
frequency  noise.  Since  the  signal  is  oversampled  at  many  times  the  Nyquist  rate,  a  digital 
low-pass  filter  may  be  used  to  remove  the  high-frequency  quantization/modulation  noise 
without  affecting  the  signal  band.  This  filtering  usually  involves  a  multi-stage 
decimation  process  since  the  output  of  the  modulator  represents  the  signal  with  the  high- 
frequency  modulation  noise  as  well  as  its  out-of-band  components  which  dominate  at  the 
lower  frequencies.  In  general,  the  smoothing  characteristics  involved  in  the  decimation 
process  require  that  the  signal  propagate  through  several  filters  and  resampling  stages. 
The  first  stage  of  decimation  lowers  the  word  rate  to  an  intermediate  frequency,  where  a 
filter  removes  the  high-frequency  modulation  noise.  A  second  low-pass  filter  is  then  used 
to  attenuate  the  out-of-band  components  before  the  signal  is  resampled  at  its  Nyquist  rate. 
As  the  signal  propagates  through  the  filters  and  resampling  stages,  the  word  length 
increases  in  order  to  preserve  the  resolution.  A  more  thorough  discussion  of  multi-stage 
decimation  and  filtering  can  be  found  in  [Ref.  1], 

The  transmission  of  coherent  light  through  optical  waveguides  has  been  of  great 
interest  ever  since  the  late  1960s.  Through  this  interest  emerged  the  concept  of 
integrated  optics,  in  which  wires  and  radio  links  are  replaced  by  light-waveguiding 
optical  fibers  and  conventional  electrical  integrated  circuits  are  replaced  by  miniaturized 
optical  integrated  circuits  (OICs).  Optical  integrated  components  offer  a  number  of 
advantages  over  their  electronic  counterparts.  These  advantages  include  large  bandwidth, 
use  of  optical  sources  capable  of  high-speed  switching  (which  is  necessary  for  high 
PRFs),  low  power  consumption,  improved  reliability  and  insensitivity  to  vibration  and 
EMI.  A  key  advantage  is  the  increased  characteristic  bandwidth  over  electronic 
components.  The  carrier  medium  is  a  lightwave  rather  than  an  electrical  current.  Thus 
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the  frequency  limiting  effects  of  capacitance  and  inductance  can  be  avoided.  Since 
electronic  converters  based  on  ZA  modulation  require  oversampling,  their  applicability  is 
mainly  limited  to  low  and  moderate  signal  frequencies.  For  instance,  oversampling  a  500 
MHz  wideband  RF  signal  by  a  factor  of  10  over  the  Nyquist  rate  would  require  a  sample 
rate  of  10  GHz,  which  is  difficult  to  process  in  an  all-electronic  bandlimited  system. 
Therefore  the  use  of  optical  integrated  components  provides  an  attractive  solution  to  the 
otherwise  bandlimited  electronic  ZA  architecture. 

B.  PRINCIPLE  CONTRIBUTIONS 


This  thesis  describes  a  single-bit,  integrated  optical  ZAM  approach.  The  ZA 
modulation  scheme  is  analyzed  and  reviewed  in  detail.  First-  and  second-order  single-bit, 
all-electronic  modulators  are  modeled  by  functional  blocks  and  simulated  using 
MATLAB's  SIMULINK  software  package.  In  the  integrated  optical  ZA  architecture, 
conventional  optical  components  are  utilized  to  match  each  functional  block  according  to 
its  system  transfer  characteristic.  The  parameters  of  these  optical  integrated  components 
are  then  optimized  to  achieve  results  matching  those  of  the  all-electronic  design.  The 
simulation  results  of  the  optical  integrated  ZAMs  illustrate  the  feasibility  of  the  integrated 
optical  approach. 

C.  THESIS  OUTLINE 

In  Chapter  II,  the  all-electronic  ZA  modulator  is  reviewed  with  both  first-  and 
second-order  architectures  being  discussed.  Analysis  and  simulation  results  for  both  the 
first-  and  second-order  models  are  presented.  Chapter  III  introduces  an  integrated  optical 
architecture  for  a  first-  and  second-order  ZA  modulator.  It  explains  the  optical  devices 
used  for  implementation  and  compares  them  at  a  component  level  to  the  all-electronic 


design.  Simulation  results  are  presented  and  analyzed.  Chapter  IV  details  the 
relationship  between  oversampling  and  resolution  and  discusses  in  terms  of  the 
modulator's  signal-to-noise  ratio.  Chapter  V  describes  the  computer  software  used  to 
simulate  the  ZA  modulators.  Finally,  Chapter  VI  discusses  issues  concerning  the  current 
simulation  model  and  presents  suggestions  for  future  efforts. 
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II.  ALL-ELECTRONIC,  SINGLE-BIT  IA  MODULATORS 


A.  FIRST-ORDER  ZAM 

A  sampled-data  equivalent  of  a  first-order  ZAM  is  shown  in  Figure  1.  Because 
this  is  a  sampled-data  circuit,  the  integration  is  performed  via  an  accumulator.  The 
analog  signal  is  assumed  to  be  oversampled  at  well  above  the  Nyquist  frequency.  This 
sampled  input,  xi5  is  fed  to  the  quantizer  via  the  accumulator.  The  quantized  output, 
which  can  be  modeled  as  an  approximation  of  the  quantization  error,  is  fed  back  and 
subtracted  from  the  input.  This  quantized  feedback  signal  forces  the  average  value  of  the 
quantized  output,  yj,  to  track  the  average  value  of  the  input  signal.  Any  difference 
accumulates  in  the  integrator  and  eventually  corrects  itself. 


Integrator/Accumulator 


The  quantization  error  is  subtracted  from  the  input  value  and  the  difference 
becomes  the  input  for  the  next  cycle.  After  the  process  is  repeated  many  times  at  high 
speed,  an  average  of  the  digital  outputs  occurring  in  each  sample  time  becomes  a  useful 
digital  representation  of  the  input  signal.  In  a  stable  converter,  the  oscillations  of  the 
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quantized  value  are  bounded,  that  is,  the  oscillations  have  a  limit  cycle.  In  general  this 
quantization  process  can  be  performed  over  more  than  one  quantization  level  [Ref.  1], 
By  this  process,  it  can  be  seen  that  the  speed  of  operation  obviates  the  need  for  precise 
circuit  elements.  Precision  in  the  quantization  levels  of  the  quantizer  is  not  a  stringent 
requirement  since  the  average  of  the  quantized  output,  yh  will  automatically  be  adjusted 
to  agree  with  the  sampled  input  analog  signal,  x,.  Therefore,  the  output  of  the  ZA 
modulation  process  can  provide  a  high  level  of  precision  in  the  representation  despite 
coarseness  in  the  quantization  levels. 

The  input/output  transfer  characteristics  of  the  first-order  ZAM  is  plotted  in 
Figure  2a.  The  signal  oscillates  between  the  quantized  levels  in  such  a  manner  that  its 
local  average  equals  the  average  input.  For  this  example  the  input  signal  is  ramped  with 
200  samples  with  a  +1  volt  range.  The  comparator  output  voltage  is  ±1  volt  with  the 
threshold  voltage  set  at  zero  volts.  Figure  2b  shows  the  undesirable  limit  cycles  at  the 
output  of  the  accumulator.  These  simulation  results  are  in  agreement  with  previously 
reported  predictions  for  first-order  ZA  modulators  [Ref.  1]. 


Samples 


Figure  2.  lst-order  all-electronic  ZAM. 

(a)  Plot  of  comparator  output  and  sampled  input, 
(b)  Plot  of  output  of  accumulator  stage. 
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B.  SECOND-ORDER  XAM 


Although  the  first-order  model  is  the  simplest,  the  quantization  noise  is  highly 
correlated  to  the  input,  resulting  in  excessive  limit  cycles.  Extending  the  architecture  to  a 
second-order  modulator  eliminates  a  number  of  instabilities  and  increases  the  reliability 
of  the  circuit.  However,  higher-order  designs  (greater  than  2)  suffer  from  instability  due 
to  the  undesirable  limit  cycles  (bounded  oscillations)  which  result  in  the  accumulation  of 
large  signals  in  the  integrators  [Ref.  4]. 

A  sampled-data  equivalent  circuit  diagram  of  a  second-order,  electronic  ZAM  is 
shown  in  Figure  3. 


First 

Accumulator 


Second 

Accumulator 


■»  I - 

[9>-(-F)-H  Delay 


L<gj-|  Delay 


t 


-4 


Figure  3.  Block  diagram  of  a  2nd-order  all-electronic  ZAM. 


The  first  accumulator,  which  embeds  the  delay  in  the  feed-backward  path,  has  a  transfer 
function  given  by 


H,(z)  = 


A 


1  -  Bz~ 


0) 


where  the  coefficients  A  and  B  are  the  gains  of  the  system.  The  second  accumulator  stage 
embeds  the  delay  in  the  feed-forward  path.  Its  transfer  function  is  given  by 


(2) 


H2{z) 


Cz~] 


1  -  Dz 


where  coefficients  C  and  D  are  the  loop  system  gains.  For  this  example  all  coefficient 
values  are  ideally  set  at  unity.  As  in  the  1  st-order  ZAM,  the  comparator  output  voltage  is 
±1  volt  with  the  threshold  voltage  set  at  zero  volts. 

The  response  of  the  second-order  ZAM  is  illustrated  in  Figure  4.  The  input  used 
for  simulation  consists  of  200  data  samples  ramped  from  -1  volt  to  +1  volt  in  increments 
of  0.01  volts.  The  duty  cycles  of  the  quantizer  output  are  weighted  toward  the  average 
value  of  the  input.  That  is,  at  the  start  of  the  ramped  input,  the  duty  cycles  are  weighted 
toward  the  bottom-level  quantization  (Figure  4a).  Towards  the  center  of  the  input,  the 
duty  cycles  are  at  about  50%.  At  the  high  end  of  the  ramp,  they  are  weighted  toward  the 
top-level  quantization.  Figure  4b  shows  the  signal  value  at  the  output  of  the  accumulator 
stages.  From  the  output  of  accumulator  #1,  it  can  be  easily  seen  that  the  output  is 
oscillating  about  the  ramped  input  range  of  -1  to  +1  volt.  Results  of  the  second-order 
modulator  illustrate  how  a  second  feedback  loop  attenuates  the  excessive  limit  cycles 
(due  to  high  correlation  of  the  quantization  noise)  found  in  the  first-order  modulator  [Ref 
!]• 

Although  not  part  of  the  current  modulator  investigation,  a  low-pass  filter  can 
then  be  used  to  resample  the  quantized  signal  at  the  Nyquist  rate.  This  serves  to  eliminate 
any  out-of-band  quantization  noise.  It  also  determines  the  ratio  of  the  sampled,  quantized 
outputs  over  the  Nyquist  interval.  This  average  value  proves  to  be  highly  representative 
of  the  input  value. 

This  chapter  provided  an  overview  of  current  sigma-delta  architecture.  The 
design  and  simulation  of  first-  and  second-order,  all-electronic  modulators  were  analyzed 
and  discussed.  These  results  provide  a  baseline  for  comparison  in  the  design  and 
simulation  of  an  integrated  optical  sigma-delta  modulator,  which  is  discussed  in  the 
following  chapter. 
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Samples 

(b) 


Figure  4.  2nd-order  all-electronic  ZAM. 

(a)  Plot  of  comparator  output  and  sampled  input, 
(b)  Plot  of  input  at  first  and  second  accumulator  stages 
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III.  INTEGRATED  OPTICAL,  SINGLE-BIT  XA  MODULATORS 

A.  FIRST-ORDER  XAM 

A  block  diagram  of  a  first-order,  integrated  optical  ZAM  is  shown  in  Figure  5.  In 
applying  optical  integrated  components  to  a  ZA  architecture,  a  first-order  model  is  first 
simulated.  In  the  integrated  optical  design,  laser  pulses  from  a  mode-locked  laser  are 
used  to  oversample  the  RF  signal.  Mode-locked  lasers  are  capable  of  providing  a  high 
pulse-repetition-frequency,  narrow  pulsewidths  and  jitter  times  on  the  order  of  200  fs.  In 
order  to  gain  a  better  understanding  of  the  model,  the  integrated  optical  components  used 
are  described  in  detail. 


RF  Signal 


Figure  5.  Block  diagram  of  a  lst-order,  integrated  optical  ZAM. 


1.  Mach-Zehnder  Interferometer 

The  Mach-Zehnder  interferometer  (MZI)  is  used  to  efficiently  couple  the 
wideband  RF  signal  into  the  optical  domain.  It  also  serves  to  subtract  the  feedback  signal 
from  the  input  signal.  Figure  6  shows  a  schematic  diagram  of  a  MZI. 


+  V 


Figure  6.  Schematic  diagram  of  a  Mach-Zehnder  interferometer  in  a  push-pull 

configuration. 

The  input  pulse  is  split  into  equal  components,  each  of  which  propagates  over  one 
arm  of  the  interferometer.  The  optical  paths  of  the  two  arms  are  equal.  If  no  phase  shift 
is  introduced  between  the  interferometer  arms,  the  two  components  combine  in  phase  at 
the  output  and  continue  to  propagate  undiminished.  For  the  current  design,  a  three- 
electrode  configuration  is  used  to  achieve  a  push-pull  phase  change  [Ref.  5],  The  push- 
pull  effect  increases  the  phase-change  efficiency  of  the  device.  This  configuration  is 
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utilized  here  to  also  subtract  the  feedback  signal  from  the  next  input  value.  In  order  to 
take  advantage  of  this  push-pull  configuration,  the  feedback  voltage  polarity  from  the 
comparator  must  be  reversed.  The  transfer  function  of  the  MZI  [Ref.  6]  can  be  expressed 
as 


where 


1 nut  ^ in 


|  +  ^cos(A<j)(F)+0) 


2nn3rTLV 
A<t>(F)  =  •GXI- 


(3) 

(4) 


is  the  voltage-dependent  phase  shift  and  is  a  function  of  the  effective  index  of  the  optical 
guide  ne,  the  pertinent  electro-optic  coefficient  r,  the  interelectrode  gap  G,  the  electrical- 
optical  overlap  parameter  T,  the  DC  bias  0,  and  the  free-space  optical  wavelength  X.  The 
modulation  voltage,  V  =VRF  -  VFB,  serves  to  subtract  the  feedback  signal  from  the  next 
input  value.  is  the  next  sampled  input  voltage  and  VFB  is  the  quantized  feedback 
voltage. 

The  method  of  accumulation  involves  the  magnitude  of  the  signal  to  be 
accumulated  and  the  direction  of  accumulation.  In  the  case  of  the  first-order  ZAM,  two 
interferometers  are  used  for  the  accumulator  stage.  One  interferometer  provides  the 
magnitude  for  the  accumulator.  The  other  interferometer  is  used  to  determine  the 
direction  of  accumulation.  Figure  7  plots  the  transfer  functions  for  both  interferometers. 
Both  MZIs  map  the  input  voltage  to  a  normalized  output  intensity  between  zero  and  one 
(light  intensity  can  not  be  negative).  The  transfer  functions  are  the  same  except  for  the 
DC  bias  0,  which  is  added  to  the  phase  shift.  For  the  MZI  controlling  the  magnitude  of 
the  signal,  0  =  n.  The  MZI  controlling  the  direction  of  accumulation  has  0  =  -n/2. 

As  can  be  seen  from  the  transfer  functions,  the  output  values  for  magnitude  range 
from  0  to  0.5  and  are  symmetric  about  the  input  value  of  zero.  The  output  values  for  the 
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direction  range  from  0  to  1 .  The  accumulator  comparator  threshold  voltage  is  normalized 
at  0.5  volts.  The  detected  direction  intensity  from  the  MZI  is  compared  to  the  normalized 
threshold  to  determine  whether  the  intensity  from  the  magnitude  MZI  accumulates 
upward  or  downward.  The  recirculating  fiber  lattice  structure  accumulates  downward  if 
the  output  of  the  interferometer  is  less  than  0.5  and  upward  for  values  greater  than  0.5. 
Thus  the  detector,  comparator  and  optical  recirculator  serve  to  function  as  an 
accumulator. 


Figure  7.  Transfer  characteristic  of  the  Mach-Zehnder  interferometer. 


2.  Fiber  Lattice  Structures 


Fiber-optic  lattice  structures  incorporating  single-mode  fibers  and  directional 
couplers  are  used  to  instrument  the  accumulators.  These  fiber  structures  can  be  used  to 
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perform  various  frequency-domain  functions  such  as  matrix  operations  and  frequency 
filtering  [Ref.  7].  Two  structures  basic  to  fiber  signal  processing  include  the  two-coupler 
nonrecirculating  and  the  two-coupler  recirculating  delay  line.  For  design  of  the  optical 
XA  modulator,  a  recirculating  feed-backward  fiber  lattice  structure  is  utilized  for  the 
accumulator.  The  generalized  form  of  a  two-coupler,  4-port  recirculating  fiber  delay  line 
is  shown  in  Figure  8.  Moslehi  et  al.  described  the  z-transform  transfer  matrix  of  this 

Adjustable  Directional  Coupler 


XI  /  Y2 


Figure  8.  Block  diagram  of  a  recirculating  fiber  optic  lattice  structure. 


recursive  structure  [Ref.  7].  The  inputs  and  outputs  are  related  to  each  other  by 


'Yx{z)\(Hn{z)  Hn{z)  YX,(z)> 
J2(z)j  [h2](z)  H22(z)j{x2(z)j 


The  transfer  functions  within  the  transfer  matrix  are 


Hu(z)  = 


a]  -(1-2(3,  )a0I|Z  1 
1  -  afl0L\Z-] 


(6) 
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H2](z)  = 


1 


1  -a0a,L,z~ 


(1  -  2 a0  -  2a,  +  4 a0a,  +  a02a 2  +  a02  -  2a0 Aa,  +  a,'  -  2 a,~a0  )L,z~ 


Hn(z)  = 


1  -aaa,L,z~ 


-a0-a,(\-2a0)L,z- 


1  -a,a0L,z 


-i 


(7) 

(8) 
(9) 


and  describe  the  use  of  the  general  structure  where  Hmn(z)  is  the  transfer  function  from 
input  X„  to  output  Ym.  Thus  H2i(z)  relates  the  X,  input  and  Y2  output  and  HI2(z)  relates 
the  X2  input  and  Y,  output.  The  parameters  a0  and  a/  are  the  intensity  coupling 
coefficients  of  the  directional  couplers  and  L,  is  the  loop  intensity  transmittance  of  the 
system.  For  convenience  Lt  is  assumed  to  be  one  (no  losses  in  the  system).  The 
accumulator  comparator  voltage  is  then  used  to  bias  the  directional  coupler  (a/)  in  order 
to  create  a  phase  change  between  the  two  pulses  in  order  to  perform  the  accumulation. 

The  accumulator  stage  in  the  first-order  model  has  a  transfer  function  given  by 
Equation  2.  From  the  transfer  matrix  of  the  recirculating  fiber  lattice  structure,  H,2(z) 
matches  this  form  where 


C  =  1  -  2an 


■  2a,  +  4  a0a,  +  aja,2  +  a02 


■2a, 


a,+a, 


-2ara, 


I  “0 


(10) 


and 


D  =  a0a, 


(11) 


are  the  corresponding  gain  values.  The  specific  lattice  configuration  is  shown  in  Figure 
9.  The  coupling  coefficients  a0  and  aI  represent  the  percentage  of  light  intensity  coupled 
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Figure  9.  Block  diagram  of  specific  fiber  lattice  configuration  used  for  accumulator 

stage  with  transfer  function  HI2(z). 

and  therefore  is  bounded  between  0  and  1.  The  desired  values  for  C  and  D  would  be 
unity.  However,  the  two  equations  work  against  each  other  simultaneously  as  shown  in 
Figure  10.  Here  a0  =  0.3  and  a,  varies  from  0  to  1 .  Values  near  the  intersection  of  C  and 


Figure  10.  Plot  of  accumulator  gains  C  and  D  as  functions  of  a0  and  aj. 
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D  provide  the  best  results  {at  =  0.5,  C  =  0.122,  D  =  0.15).  In  order  to  compensate  for  the 
small  value  of  C,  an  optical  amplifier  with  a  gain  of  100  is  placed  just  prior  to  the  fiber 
lattice  structure  in  the  accumulator  stage  (see  Figure  5).  Figure  11  plots  the  transfer 
characteristic  of  the  first-order,  integrated  optical  ZAM.  It  is  apparent  that  the  first-order 
system  is  not  stable  enough  for  accurate  conversion  of  the  signal  due  to  excessive  limit 
cycles. 


Samples 


Figure  11.  Transfer  characteristic  of  1  st-order  electro-optic  ZAM. 
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B.  SECOND-ORDER  XAM 


The  block  diagram  for  a  second-order,  integrated  optical  LAM  is  shown  in  Figure 

12. 


Figure  12.  Block  diagram  of  a  2nd-order  electro-optic  LAM. 


The  first  stage  (Accumulator#  1)  has  the  transfer  function  described  by  Equation  1  and 
uses  H2,(z)  given  by  Equation  7.  The  specific  fiber  lattice  structure  configuration  used 
for  the  first  accumulator  stage  is  shown  in  Figure  13. 
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Figure  13.  Block  diagram  of  specific  fiber  lattice  configuration  used  for  accumulator 

stage  with  transfer  function  H2j(z). 

The  coupler  coefficients  are  set  ideally  at  unity,  thus  the  gains  A  and  B  are  also 
set  at  unity.  The  second  stage  (Accumulator#2)  is  identical  to  the  accumulator  stage  in 
the  first-order  XAM  (see  Figure  9).  The  coupler  coefficients  are  again  a0  =  0.3  and  a,  = 
0.5;  however,  the  gain  of  the  optical  amplifier  is  now  15.  The  values  for  a0  and  a,  were 
found  to  best  optimize  the  accumulator  gains  C  and  D.  (Appendix  A  provides  simulation 
results  for  several  other  values  for  a0  and  a,.) 

Simulation  results  for  the  second-order,  integrated  optical  XAM  are  plotted  in 
Figure  14.  The  average  value  of  the  quantizer  output  can  be  seen  to  track  the  average 
value  of  the  ramped  input  as  shown  in  Figure  14a.  The  output  of  the  interferometer  can 
be  seen  to  oscillate  about  the  ramped  input.  Figure  14b  plots  the  intermediate  signal 
values  at  the  input  of  the  MZIs  in  the  accumulator  stages. 

This  chapter  detailed  the  design  and  simulation  for  an  integrated  optical  XAM. 
Results  for  both  first-  and  second-order  modulators  were  analyzed.  These  results 
compared  favorably  to  those  of  the  all-electronic  design  and  demonstrated  the  feasibility 
of  the  integrated  optical  approach.  The  next  chapter  discusses  the  effects  of  signal-to- 
noise  ratio  and  oversampling  ratio  of  on  the  bit  resolution  of  XAMs. 
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Figure  14.  2nd-order  electro-optic  ZAM. 

(a)  Plot  of  comparator  output  and  sampled  input, 
(b)  Plot  of  input  at  first  and  second  accumulator  stages 


IV.  OVERSAMPLING  RATIO  AND  RESOLUTION 


The  quantization  utilized  in  LAM  introduces  noise  in  the  modulator.  The 
quantization  error  e  is  treated  as  white  noise  having  probability  of  lying  anywhere  in  the 
range  ±A/2,  where  A  is  the  level  spacing  (normalized  units)  between  quantized  levels. 
Its  mean  square  value  is  given  by  [Ref.  1] 


(12) 


The  oversampling  ratio  (OSR),  defined  as  the  ratio  of  the  sampling  frequency  f s  to  the 
Nyquist  rate  2/0 ,  is  given  by  the  integer 


OSR  = 


2/0 


(13) 


The  noise  power  in  the  signal  band  can  be  shown  to  be  [Ref.  1  ] 

K  04) 

It  is  evident  that  oversampling  reduces  the  in-band  rms  quantization  noise  n0  by  the 
square  root  of  the  OSR. 

The  feedback  loops  in  the  ZAM  help  shape  the  spectrum  of  the  modulation  noise 
by  moving  most  of  the  noise  outside  the  signal  band.  The  filters  used  in  the  loops  reduce 
the  net  noise  in  the  signal  band.  The  ZAM  subtracts  the  previous  value  of  the 
quantization  error  from  the  present  error.  In  the  case  of  two  feedback  loops,  the 
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modulation  noise  becomes  the  second  difference  of  the  quantization  error.  The  signal-to- 
noise  ratio  for  a  second-order  ZAM  can  be  predicted  from  [Ref.  1] 


^  =  ^(°SRT*.  (“> 

For  a  second-order  ZAM,  the  signal-to-noise  ratio  (SNR)  increases  at  15  dB/octave  and  6 
dB/bit.  Thus  for  OSR  =128  and  A  =  2,  n0  =  -97  dB.  For  a  signal  strength  of  one  (i.e.,  0 
dB),  SNR  =  -n0  =  97  dB.  At  6  dB/bit,  this  translates  to  16  bits  of  resolution. 

For  example,  the  resolution  of  the  second-order,  integrated  optical  ZAM  can  be 
determined  using  the  above  method.  Assuming  the  200  input  data  points  are  sampled 
over  one  second  (f0  =  1  Hz),  the  oversampling  ratio  then  becomes 


OSR  = 


200 
2(1  Hz) 


=  100. 


(16) 


Normalizing  to  a  signal  strength  of  one,  SNR  =  -w0  =  91  dB,  which  translates  to  15  bits 
of  resolution. 

When  discussing  ADC  architectures,  it  is  important  to  analyze  how  bit  resolution 
is  determined.  This  chapter  discussed  in  detail  how  bit  resolution  of  ZAMs  is  affected  by 
the  SNR  and  oversampling  ratio.  Bit  resolution  is  proven  to  be  a  function  the 
oversampling  ratio.  The  following  chapter  describes  the  computer  software  utilized  to 
simulate  modulator  designs. 
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V.  SIMULATION  TOOLS 


In  simulating  the  ZA  modulators,  MATLAB's  SIMULINK  toolbox  was  used. 
SIMULINK  is  a  program  for  simulating  dynamic  systems.  It  is  an  extension  of 
MATLAB  and  uses  block  diagrams  to  represent  dynamic  systems.  The  purpose  of  this 
chapter  is  to  familiarize  the  reader  with  SIMULINK  and  how  simulation  results  were 
obtained.  Figure  15  shows  the  system  block  diagram  for  the  first-order,  all-electronic 
ZAM.  The  simulation  variables  input ,  suml,  acl  and  output  are  are  used  for  plotting 
simulation  results.  Appendix  B  lists  the  MATLAB  functions  and  variables  used  in 
SIMULINK. 


Figure  15.  SIMULINK  diagram  of  lst-order,  all-electronic  ZAM  (filename  ell.m). 

The  ramped  input  is  a  MATLAB  vector  (elec  in.mat)  of  200  data  values  from  -1  to  +1 
incremented  at  0.01.  This  input  is  used  for  all  simulations.  The  unit  delay  is  described 
by  the  z-domain  transfer  function  1/z  and  is  embedded  in  a  feedback  loop.  The  relay 
block  is  a  binary  switch  that  outputs  a  +1  if  the  input  is  above  zero  and  -1  if  the  input  is 
below  zero. 

The  second-order,  all-electronic  ZAM  is  modeled  in  much  the  same  fashion, 
except  that  a  second  feedback  loop  is  added  (see  Figure  16).  The  system  gains  are  set  at 
unity  for  simulation  purposes. 
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Figure  16.  SIMULINK  diagram  of  2nd-order,  all-electronic  ZAM  (filename  el2.m). 


The  SIMULINK  block  diagram  used  to  model  a  first-order,  optical  integrated 
ZAM  is  shown  in  Figure  17.  The  input  vector  (optic  in.  mat)  is  identical  to  that  of  the  all- 
electronic  model. 


Figure  17.  SIMULINK  diagram  of  lst-order,  optical  integrated  ZAM  (filename  opl.m). 

The  MZI  block,  unmasked  in  Figure  18,  is  used  to  subtract  the  two  inputs  and  map  the  the 
input  according  to  the  interferometer  transfer  characteristic  (Equation  3).  The  MATLAB 
code  for  the  function  describing  the  interferometer  (MZI)  is  listed  in  Appendix  B.  After 
the  signal  is  mapped  to  a  value  between  zero  and  one,  another  MATLAB  function  (v-0.5) 
is  used  to  subtract  0.5  from  the  input  (v).  The  subtraction  is  done  to  normalize  the  signal 
about  zero  in  order  for  the  discrete  integrators  to  function  properly. 


26 


in_2 

Figure  18.  Unmasked  MZI  block  from  Figure  17. 

The  signal  then  passes  through  a  gain  stage  of  100  before  entering  the 
integrator/accumulator.  The  integrator  is  basically  a  filter  which  can  be  describe  by  z- 
domain  characteristics.  The  coefficients  for  this  filter  was  described  in  detail  in  Chapter 
IV.  Finally,  the  relay  used  is  the  same  switch  used  in  the  all-electronic  simulation  model. 

Figure  19  shows  the  SIMULINK  block  diagram  of  the  second-order,  optical 
integrated  ZAM.  All  functional  blocks  are  identical  to  the  first-order  model  with  the 
exception  of  the  filters  (Accum#l  and  Accum#2).  The  coefficients  for  these  discrete-time 
filters  are  described  in  Chapter  IV. 
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Figure  19.  SIMULINK  diagram  of  2nd-order,  optical  integrated  ZAM  (filename  op2.m). 

This  chapter  described  the  use  of  computer  simulation  software  (MATLAB  and 
SIMULINK)  used  to  design  and  simulate  all  ZA  modulators  in  the  thesis.  The  figures 
provide  system  block  diagrams  and  filenames  used  in  obtaining  simulation  results.  The 
final  chapter  discusses  the  current  design  and  suggest  recommendations  for  future  efforts. 
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VI.  CONCLUSIONS 


The  LA  oversampling  A/D  modulator  architecture  uses  limit  cycles  in  quantized 
feedback  loops  to  provide  an  accurate  digital  representation  of  the  input  signal.  The 
second-order  LAM  provides  a  stable  and  robust  design  that  is  highly  tolerant  of  circuit 
imperfections  and  component  mismatch.  The  major  limitations  of  this  method  are  fast 
cycle  times  and  bandwidth.  The  use  of  fiber  optic  technology  has  the  potential  of 
eliminating  these  limitations.  An  integrated  optical  second-order  LA  architecture  allows 
the  processing  of  wideband  RF  signals.  The  integrated  optical  LAM  design  presented  in 
this  paper  is  a  fairly  straight-forward  extension  of  the  electronic  design  using  standard 
integrated  optical  devices.  Current  simulation  results  confirm  design  feasibility. 

Future  efforts  include  further  optimization  of  the  current  integrated  optical  design. 
Modifications  may  include  the  possibility  of  optimizing  the  magnitude  and  direction  of 
each  accumulation  stage  using  only  one  interferometer.  The  accumulation  stages  of  the 
integrated  optical  LAM  are  modeled  by  fiber  lattice  structures  with  similar  system 
transfer  functions.  More  rigorous  modeling  of  these  delay  line  structures  and  analysis  of 
the  optical  amplifiers  are  needed  for  future  hardware  implementation.  Since  the  output  of 
the  modulator  represents  the  input  signal  together  with  modulation  noise,  there  is  still  a 
need  to  decimate  the  modulated  signal.  A  multi-stage  decimation  is  needed  to  lower  the 
word  rate  and  remove  high-frequency  modulation  noise  before  the  signal  is  resampled  at 
the  Nyquist  rate.  Design  issues  to  be  studied  include  nonlinearities  associated  with 
interferometers,  stability  of  the  accumulator  stages,  effects  of  net  gains  in  the  feedback 
loops,  and  the  effects  of  modulation  noise  and  oversampling  frequency  (OSR)  on  bit 
resolution. 
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APPENDIX  A:  SIMULATION  RESULTS 


The  following  figures  provide  further  simulation  results  for  the  second-order, 
optical  integrated  ZAM  using  several  other  values  for  the  directional  coupler  coefficients 
(a0  and  a,).  These  results  are  provided  to  illustrate  the  sensitivity  of  the  system  to  the 
coefficient  settings. 
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Samples 


Figure  21.  Plot  of  output  and  input  with  gains  C  =  0. 1 6  and  D  =  0. 1 0. 


Samples 

Figure  22.  Plot  of  output  and  input  with  gains  C  =  0.058  and  D  =  0.14 


Volts 
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APPENDIX  B:  MATLAB  CODE  USED  IN  SIMULINK 


FILENAME  interfer.m 

%  INTERFEROMETER  TRANSFER  FUNCTION  (Equation  3) 
function  I_out=interfer(v) 


G=3.0e-6; 

wl=0.9e-6; 

n=2.205; 

r=30.8e-12; 

GAMMA=0.5; 

L=0.0020; 

V_rf=0; 

Vs=0; 


%  Gap  (m) 

%  wavelength  (m) 

%  index 

%  electro-optic  coefficient  (v/m) 
%  overlap  integral 
%  length  (m) 


phase=-pi/2; 

K=(2  *pi*nA3*r*  GAMMA  *  L)/(G  *  wl); 
I_out=(.5  +  ,5*cos(K*v  +  phase)); 


FILENAME  el_ramp.m 

%  Input  ramp  vector  for  both  electronic  and  optical  models 
limit=l; 

V_rf=[-limit:.01 :  limit]; 
t=[0:limit*200]; 

V_in=[t;V_rf]; 
input=V_in; 
save  elec  in  input 


FILENAME  ell_plot.m 

%  Plots  ell.m  output  (First-Order,  All-Electronic  Model) 

save  ell  out  suml  input  output  %  Saves  simulation  variables  for  plot 

clear 

limit=l; 
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figure(l)  %  First  figure  plots  input  and  output 
clg 

load  ell_out 

subplot(2, 1 , 1  ),plot(input,'c'),hold,stairs(output) 

set(get(gca,'Chil'), 'color', 'red') 

ylabel(' Volts') 

axis([0  limit*200  -1.5  1.5]) 

set(gca,'xtick',[]); 

subplot(2, 1 ,2),plot(suml  ,’c') 

axis([0  limit*200  -limit-1  (limit+1)]) 

ylabel('Volts'),xlabel('Samples') 


FILENAME  el2_plot.m 

%  Plots  el2.m  outputs  (Second-Order,  All-Electronic  Model) 

save  el2_out  suml  acl  sum2  ac2  input  output  %  Saves  simulation  variables  for  plot 

clear 

limit=l; 

figure(  1 )  %  F irst  figure  plots  input  and  output 

clg 

load  el2_out 

plot(input,’c'),hold,stairs(output) 
set(get(gca,'Chil'), 'color', 'red') 
ylabel(’Volts'),xlabel('Samples') 
axis([0  limit*200  -1.5  1.5]) 

figure(2)  %  Second  figure  plots  intermediate  stages 
subplot(2, 1 , 1  ),plot(sum  1  ,’r'),title(' Accumulator#  1 ') 
ylabel(' Volts') 
set(gca,'xtick',[]); 

subplot(2, 1 ,2),plot(sum2,'r'),title('Accumulator#2') 

ylabel(' Volts') 

xlabel('Samples') 
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FILENAME  opl_plot.m 

%  Plots  op l.m  outputs  (First-Order,  Optical  Integrated  Model) 

%  Save  simulation  variables  for  plot 
save  opl_out  modi  inti  sig_in  sig_out 

clear 

limit=l; 

figure(l) 

clg 

load  opl_out 

subplot(2,l,l),plot(sig_in,'c’),hold,stairs(sig_out) 
ylabel(' Volts') 
set(get(gca,'Chil'),'colorVr') 
set(gca,'xtick',[]); 

axis([0  limit*200  -limit-.5  limit+.5]) 
subplot(2, 1 ,2),plot(mod  1  ,'c') 
ylabel('Volts'),xlabel('Samples') 


FILENAME  op2_plot.m 

%  Plots  op2.m  outputs  (Second-Order,  Optical  Integrated  Model) 

save  op2_out  modi  inti  mod2  int2  sig  in  sig_out 

clear 

limit=l; 

figure(l) 

clg 

load  op2_out 

plot(sig_in,'c'),hold,stairs(sig_out) 
set(get(gca,'Chil'),'color','red') 
ylabel('Volts'),xlabel('Samples') 
axis([0  limit*200  -limit-. 5  limit+.5]) 
figure(2) 

subplot(2,l,l),plot(modl,'r'),title('(Magnitude)  MZI  #1') 

ylabel(' Volts') 

set(gca,'xtick',[]); 

subplot(2,l,2),plot(mod2,'r'),title('(Magnitude)  MZI  #2') 
ylabel(' Volts') 
xlabel('  Samples') 
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FILENAME  mzi.m 

%  Plots  MZI  transfer  characteristics 

clear 


G=3.0e-6; 

wl=0.9e-6; 

n=2.205; 

r=30.8e-12; 

GAMMA=0.5; 

L=0.0040; 


%  Gap  (meter) 

%  wavelength  (meter) 

%  index 

%  electro-optic  coefficient  (v/m) 
%  overlap  integral 
%  length  (m) 


V_rf=0; 

Vs=0; 

sign=-pi/2; 

mag=pi; 


vpi=G*wl/(2*L*nA3*r*GAMMA) 


v=[-1.5:.05:1.5]; 


K=(2*pi*nA3*r*GAMMA*L)/(G*wl); 
for  i=  1:61 

I_sign(i)=(.5  +  .5*cos(K.*v(i)  +  sign)); 
end 

for  j=l:61 

I_mag(j)=(.5  +  ,5*cos(K.*v(j)  +  mag)); 
end 

clg 

plot(v,I_sign,'r+'),hold,plot(v,I_mag,'co') 
xlabel('Normalized  Input  Signal  (Volts)') 
ylabel('Normalized  Output  Intensity') 
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FILENAME  Iattice.m  (Equations  10  and  11) 

%  Plots  coefficients  for  Lattice  transfer  function 

clear 

inc  =  0.01; 

C  =  zeros(  1  ,inc*  1 0000); 

D  =  zeros(  l,inc*  10000); 

aO  =  .3; 

LI  =  1; 

for  al  =  .01:inc:1.00 

C(al  *  1 00)  =  ( 1  -2*a0-2*al+4*al  *a0+al  A2*aOA2+aOA2-2*aOA2*al ... 
+al  A2-2*al  A2*aO)*L  1 ;  %  These  two  lines  constitute  one  equation 

D(al*100)  =  a0*al*Ll; 
end 

clg 

plot(C,'c-'),hold,plot(D,'r-’) 
xlabel('Coupler  Coefficient  al') 
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